Effects of morphology on phonons of nanoscopic silver grains 
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The morphology of nanoscopic Ag grains significantly affects the phonons. Atomistic simulations 
show that realistic nanograin models display complex vibrational properties. (1) Single-crystalline 
grains. Nearly-pure torsional and radial phonons appear at low frequencies. For low-energy, faceted 
models, the breathing mode and acoustic gap (lowest frequency) are about 10% lower than predicted 
by elasticity theory (ET) for a continuum sphere of the same volume. The sharp edges and the 
atomic lattice split the ET-acoustic-gap quintet into a doublet and triplet. The surface protrusions 
associated with nearly spherical, high-energy models produce a smaller acoustic gap and a higher 
vibrational density of states (DOS) at frequencies u < 2THz. (2) Twined icosahedra. In contrast 
to the single-crystal case, the inherent strain produce a larger acoustic gap, while the core atoms 
yield a DOS tail extending beyond the highest frequency of single-crystalline grains. (3) Mark's 
decahedra, in contrast to (1) and (2), do not have a breathing mode; although twined and strained, 
do not exhibit a high-frequency tail in the DOS. (4) Irregular nanograins. Grain boundaries and 
surface disorder yield non-degenerate phonon frequencies, and significantly smaller acoustic gap. 
Only these nanograins exhibit a low-frequency v 2 DOS in the interval l-2THz. 



I. INTRODUCTION 

Recent experiments have probed the size dependence of 
the electronic-energy relaxation rate in nanoscopic metal 
grains Perhaps the most striking aspect of these ex- 
periments is the finding of Arbouet and coworkers which 
shows that the electron-phonon relaxation rate r ep does 
not depend on the environment of the grains X Gold and 
silver grains suspended in solution present the same size 
dependence of r ep than those embedded in a solid matrix 
(e.g., A1 2 3 ). 

To explain the data, current models for the energy re- 
laxation of hot electrons due to electron-phonon coupling 
in nanoscopic grains rely on phenomenological parame- 
ters and simplifications. (For example, replacing the vi- 
brational density of states with the Debye model.) One 
has to keep in mind, however, that the experiments in- 
volved several materials (Ag, Au, Pt) with a wide range 
of grain sizes, where a myriad of stable low-energy mor- 
phologies has been observed and proposed Thus, to 
understand the electron-phonon coupling in nanoscopic 
metal grains it would be desirable to have a description 
of their vibrational properties that includes the grain's 
morphology: facets, grain boundaries, and surface disor- 
der and defects. 

There are papers in the literature that have addressed 
the vibrational density of states of nanoscopic metal 
grains using ideal, spherical models&2ii£ Unfortunately, 
the use of these idealized models limits a detailed un- 
derstanding of the vibrational properties and the role 
of morphology. Furthermore, the idealized models ne- 
glect the influence of grain boundaries, and surface dis- 
order and defects on the vibrational properties of the 
experimental nanograins. 11 In this work we focus on re- 
alistic, nanoscopic Ag grains containing 150-1400 atoms 



and atomistically simulate their vibrational properties, 
by calculating the phonon frequencies and displacement 
vectors. We show that morphology introduces a high 
degree of complexity in the phonon spectra; total and 
partial vibrational density of states; and phonon local- 
ization. 

Our predictions, which should be general for 
nanoscopic metals, show the following prominent fea- 
tures, (a) Low-energy, single-crystalline grains present 
nearly-pure torsional and radial phonon modes. The fre- 
quencies of the breathing mode and the acoustic gap 
(lowest frequency) are nearly 10% lower than predicted 
by elasticity theory (ET) for a continuum sphere .iSiiiii^ 
The sharp edges and atomic lattice of the grains lead 
to the splitting of the acoustic gap quintet predicted by 
ET in a doublet and triplet, where the magnitude of the 
splitting depends on the relative number of atoms at the 
boundary of different facets. High-energy, ideal spher- 
ical models present regular protrusions on the surface. 
When compared to a faceted grain of the same size, these 
protrusions lead to a smaller acoustic gap, which is not 
necessarily degenerate; and to a higher total vibrational 
density of states (DOS) at frequencies v < 2 THz. 

(b) As in the case of single-crystalline nanograins, 
twined icosahedra have a breathing mode. Strain in 
these grains leads to a higher acoustic gap, and a high- 
frequency tail in the DOS that extends beyond the high- 
est frequency in single-crystalline grains. On the other 
hand, twined and strained Mark's decahedra do not have 
a breathing mode, neither exhibit a high-frequency tail 
on DOS. 

(c) Irregular nanograins with grain boundaries and sur- 
face disorder do not have degenerate phonon frequen- 
cies, and the acoustic gap is significantly reduced. These 
nanograins are the only ones to exhibit v 2 DOS in the 
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interval l-2THz. The extent of this region depends on 
the nature of disorder. 

To summarize, our simulations show the complexity 
of the vibrational properties of nanoscopic metals, and 
point out similarities between the vibrational density 
of states of nanograms with grain boundaries and sur- 
face disorder and that of massive nanocrystalline sam- 
ples. This work is a required step towards understanding 
the size dependence of the electron-phonon coupling in 
nanoscopic metal grains. 



II. NANOGRAINS' MORPHOLOGY AND 
ENERGETICS 



The morphology of a nanograin is specified by the num- 
ber of atoms N, shape and atomic arrangement. The 
shape refers to the symmetry and facets, while the atomic 
arrangement pertains to the nanograin's crystalline char- 
acteristics: single-crystalline, twined, surface-disordered 
and defective. Here, we classify the model nanograins — 
with sizes that range from N= 150-1400 — in three broad 
classes: 

(I) Single-crystalline nanograins derived from the face- 
centered-cubic lattice of Ag with facets oriented in the 
[111], [001], and [101] directions (depending on size) . Fig- 
ure ^a) illustrates class I grains. TO and TO + origi- 
nate from different degrees of truncation of an octahedron 
(Oh). TO's are regular in the sense that the number of 
atoms in every edge is the same and the [111] and [001] 
facets form regular hexagons and squares, respectively. 
TO + represents a group of grains that result from re- 
moving the vertices of the Oh first and then [001] facets 
until reaching the TO. The latter protocol leads to grains 
that have more atoms in the edges connecting [111] facets 
than in the [111]/ [001] boundary. Model C results from 
truncating cuboctahedra. In these grains, [001] facets are 
larger than the corresponding facets in TO + , while the 
geometry of [111] facets is preserved. On the other hand, 
by simply removing the atoms lying in the [lll]-facets 
boundary, Model C leads to Model H. By construction, 
Model H nanograins show [101] facets. 

(II) Twined icosahedra (Ih) and its variants, 15 and 
Mark's decahedra (M-Dh). Examples of Ih and M-Dh ap- 
pear in Fig. n Both these grains are described in detail 
in the literaturei^S^i Ih's are strained and feature [111] 
facets, 12 vertices, 20 internal grains, and six 5-fold axis; 
one of them appears indicated with a "5" in Fig. ^ I n 
contrast, M-Dh's are less strained, have multiple faceting 
and a single 5-fold axis. 

(III) Complex, partially disordered multigrain 
nanograins that have grain boundaries, surface disorder, 
and defects; see Fig. QIc). Simulated annealing is used 
to generate these grains. The most common surface 
defects are vertex vacancies [shown in grey in Fig. ^c)], 
terraces, and stacking faults. In most cases, the grain 
boundaries lead to 5-fold-like axis and a high number 
(about 10) of internal grains. However, we have also 
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FIG. 1: Classes of Ag nanograins illustrated with selected 
specimens. "5" [(b)] and "5 '" [(c)] indicate, respectively, a 
5-fold and 5-fold-like axis. Number of atoms (N) is indi- 
cated, (a) Single-crystal grains with [111] and [001] facets in 
the regular (TO) and irregular (TO + ) truncated octahedron 
geometry. Compared to TO + , Model C has bigger [001] facets 
and similar [lll]-facet geometry. Model H grains present [110] 
facets, (b) Twined icosahedron (Ih) and Mark's decahedron 
(M-Dh). (c) Multigrain nanograins exhibit grain boundaries 
and surface disorder; gray atoms highlight vertex vacancies. 
Spherical nanograins result from decorating regular, faceted 
specimens. Gray atoms decorating [111] facets of TO [see (a)] 
yield the shown "sphere." 



encountered multigrains with a single twin boundary or 
internal stacking fault. 

We also include in class III the so-called spherical 
nanograins [Figure ^c)]. Despite their unrealistic mor- 
phology, these grains have been the model of choice 
in recent simulations of the vibrational properties of 
nanoscopic metal grains. 8,10 Although single-crystalline, 
these grains have small protrusions (grey atoms in the 
figure) on the [111] and, depending on size, [001] facets 
of otherwise compact specimens. These protrusions act 
as surface defects, despite of being highly regular. (In- 
deed, one may consider the compact nanograins to be 
decorated by the protrusions. These surface defects con- 
tribute to the higher energies of the spherical nanograins 
relative to other structures and are removed by simulated 
annealing starting at high temperatures. 

Figure [2] shows the effect of morphology on the to- 
tal energy £ for class I and II. Ih are the most stable 
grains up to TV ~ 600, where the cohesive energies of 
Ih, M-Dh and TO are effectively degenerate. For larger 
grains, strain dominates over efficient faceting and single- 
crystalline grains become more stable. TO + splits into 
two bands, indexed by the number of atoms in [lll]-facet 
boundaries; the lower and upper bands corresponding to 
5 and 7 atoms, respectively. The lower and upper band 
merge at N — 1289 and N = 1975, respectively, with TO 
energies. Model H tracks the TO + upper band. Mark's 
decahedra show significant fluctuations in the size depen- 
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FIG. 2: Morphology effects on the total energy £ of grains 
in Fig. 0^a) and (b). £ is measured relative to bulk (es = 
2.85 eV). Twined icosahedra (Ih) are the most stable struc- 
tures from 150 to N ~ 600; where the cohesive energies of 
Ih, M-Dh, and TO are nearly degenerate. The energy of TO + 
nanograins separates in two bands (see Sec. [2}. The upper 
band follows the progression of Model H. Increased faceting is 
responsible for the noticeable scattering in M-Dh's energy. 

TABLE I: Total energy of multigrain and spherical Ag grains 
vs number of atoms N. Here, eb — 2.85 eV is the bulk cohe- 
sive energy. 
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dence of £ due to increased faceting; which permits gen- 
crating M-Dh grains with significantly different sphericity 
even if the number of atoms is slightly changed^ 

Considering the energetics of class III, Table [I] lists 
the total energies of multigrain and spherical nanograins. 
These energies exceed by a considerable amount those 
of class I and II. Nonetheless, multigrain nanograins are 
non-equilibrium structures that are likely to exist under 
experimental conditions. In addition, the recent simu- 
lations of the structure of ultra-small Al clusters (Man- 
nincn et at, in Ref. ITsf) show that complex, multigrain 
structures are competitive with those of class I and II. 
Spherical grains, on the other hand, correspond to an 
idealization of realistic metal grains and are presented 
here only for comparison purposes. 



III. VIBRATIONAL PROPERTIES 

We show that nanograin morphology strongly af- 
fects the vibrational properties. By using the meth- 
ods described in the Appendix, we analyze the effect of 



nanograin's morphology on the phonon frequency spec- 
trum, vibrational density of states (DOS), and phonon lo- 
calization. First, we describe general features of phonon 
spectra followed by a comparison between class I and 
spherical grains with class II and multigrains. Second, we 
contrast the DOS of TO with class II and III nanograins. 
Finally, we focus on N=459 nanograins in class I and III 
to illustrate important features of the phonon localiza- 
tion. 



A. Phonon frequencies 

General features. The finite size of the grains lead to 
discrete phonon spectra. The acoustic gap, which is de- 
fined as the value of the lowest phonon frequency, de- 
creases with size due to confinement. Our simulations 
show that the acoustic gap roughly scales as iV -1 / 3 , as 
predicted by elasticity theory (ET) for a finite grain^ 
and its magnitude depends on morphology. Each grain 
in class I as well as icosahedra (Ih) have a nearly-radial, 
breathing mode. In contrast, this mode does not appear 
in Mark's decahedra (M-Dh) or multigrains (class III). 
When compared with predictions of elasticity theory for 
a solid sphere with the same volume as the nanograins, 
our simulations lead to frequencies of the breathing mode 
(i?i) and acoustic gap about 10% smaller. For R±, ET 
predicts — 2.825 (ci/2nR), where R is the radius of 
the sphere and c; = 3686 m/s is the average longitudi- 
nal sound velocity of bulk silver The ET-acoustic-gap 
equals 0.401 with 5-fold degenerac y 13 i 14 

Class I and spherical nanograms. For TO, TO + , and 
Model C in class I, and for the spherical grains in class 
III, Figure [21 shows phonon frequencies up to the first 
nearly-radial, breathing mode. Main features are the fol- 
lowing, (i) All grains show frequencies that are either 
non-degenerate or 2-fold and 3-fold degenerate, (ii) The 
five lowest frequencies split into a triplet and a doublet; 
a gap separates this cluster from other frequencies. This 
triplet-doublet cluster originates from symmetry break- 
ing of the ET-acoustic-gap quintet. The size dependence 
of the small split-off separation between the triplet and 
doublet depends on the morphology. It decreases for TO, 
and in TO + , N=459 the triplet and doublet even switch 
order. Model C shows a more remarkable size depen- 
dence: The split-off separation remains nearly constant 
and similar to that of N=201 TO grain. We attribute this 
behavior to the fixed number of atoms (3) in the [111]- 
facet boundaries, (iii) In TO + and in 586 and 1289 TO 
grains, a distinguishable, second group of frequencies ap- 
pears; a triplet precedes a doublet and two triplets follow. 
The separation between multiplets in these complexes 
also depend upon the size of the nanograin. In contrast, 
no visible gap separates these frequencies in Model C and 
201, TO. At frequencies higher than 0.5i/j? no sizeable 
gaps appear in any grain, (iv) In all grains, near to the 
lowest triplet-doublet, we find a low- frequency nearly- 
torsional edge mode T± ; morphology dramatically affects 
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FIG. 3: Single-crystalline Ag grains' scaled phonon frequen- 
cies v up to the first nearly-radial, breathing mode (Ri). 
Vp = 2.825 (ci/2-kR) (see Sec. IIII Al . Frequencies are ei- 
ther nondegenerate (short dash) or 2-fold (medium dash) and 
3-fold (long dash) degenerate frequencies are present. In TO, 
TO + and Model C, the low-lying phonon frequency structure is 
qualitatively similar except for Ti ; a nearly-torsional, edge vi- 
bration (sec Fig. T2 labels a 2-fold degenerate twist mode. 
Asterisks (*) indicate nearly degenerate doublet-triplet clus- 
ters. Surface protrusions yield different low-frequency phonon 
structure in spherical grains. 



the size dependence of its frequency (Fig. In TO, 
TVs frequency remains almost constant while in TO + it 
decreases with the surface. In Model C, it scales as 
and is locked with the lowest five frequencies. 

The ideal spherical nanograins present phonon fre- 
quencies that are nondegenerate as well 2- and 3-fold 
degenerate. Contrary to class I grains, the low-lying 
frequencies present a less clear pattern. These grains, 
N=225 and 459, may be viewed as the decorated version 
of the 201 TO (c.f. Fig. HJ and 405 TO+, respectively. 
The coupling between the vibration of the protrusions 
and the atoms in the faceted grain that they decorate 
reduces the acoustic gap. This coupling also destroys the 
similarities between the frequency structure in 201 and 
405. (See Fig. H) 

For N=586 TO nanogram, Fig. Ufa) shows represen- 
tations of the displacement vectors of 3 selected phonon 
modes. As mentioned above, T± is nearly torsional and 
the atoms surrounding the edges have the higher ampli- 
tude. T2 is one of the two nearly-torsional, 2-fold gener- 
ate twist mode. In contrast to Ti, T^'s frequency depends 
less strongly on morphology (see Fig. |3J). T\ and T2 also 
appear in spherical grains. i?i is the non-degenerate, 
nearly-radial breathing mode. It should be noted that 
TO and Ih grains have another nearly-torsional mode T3 
(not shown), which is an excited, non-degenerate twist 
mode. This mode does not appear in any other single- 
crystalline grain or Mark decahedron, as only TO's and 
Ih's have the same number of atoms at the every [111]- 




?, R 1 
mm 1 J y >m» 

inkMt > > > > J : » 
■ ■ > •» » » > ) x : : ; 

war* 1 i > ww w 

Pilfer 




FIG. 4: (a) Displacement vectors of selected phonon modes 
in the N=586, TO nanogram: (i) Ti, a non-degenerate mode 
in which the atoms surrounding the edges have a significant 
amplitude, (ii) T2, one of the doubly-degenerate twist modes, 
and (iii) Ri, the first nearly-radial, breathing mode, (b) Dis- 
tribution functions (1° bin size) of the angle between the equi- 
librium vector position of an atom — measured from the center 
of the grain — and its displacement vector. The distributions 
peak around 90° for Ti and T2, and 10 (instead zero) for Hi, 
respectively. For larger grains the distributions peak around 
2-4°. 



facet boundary, [c.f. Figure HJa) and (b)] 

For Ti, T 2 , figure |H[b) shows histograms of the an- 
gles between the displacement vector of atom i (u 1 ) and 
its equilibrium position (i? 1 ); measured from the center of 
the grain. [The angle is defined mod(7r).] The purpose is 
to show the extent to which the crystalline field and edges 
of the grains couple angular and radial degree's of free- 
dom. Note that in an ideal elastic sphere, torsional and 
breathing modes have angular- and radial-only displace- 
ments, respectively^ The values for the average angle 9 
and the standard deviation a of the histograms in Fig. 
Hare the following: 9 Tl = 89.51°, a Tl = 12.53°; 9 T2 = 
89.95°, ctt 2 = 14.46°; and Rl = 10.64°, a Rl = 5.96°. 
In general, as size increases, for nearly torsional and ra- 
dial modes, 9 approaches 90° and 0°, respectively; and a 
decreases. 

Class II and multigrain nanograins. Figure[5fa) shows 
the phonon frequencies up to the first nearly-radial, 
breathing mode of twined icosahedra (Ih). Panel (b) 
shows phonon frequencies up to Vp for Mark's decahedra 
(M-Dh) and multigrain samples. (See Fig. QJ) Prominent 
features are the following: 

(i) Strain in the icosahedra leads to clusters of non- 
degenerate frequencies. At higher frequencies a few 2- 
fold degeneracies appear. The acoustic gap is higher 
than in class I and spherical grains. It also has a triplet- 
doublet structure with higher split-off and no degenera- 
cies, (ii) Mark's decahedra, like Ih's, only show a few 
doubly degenerate frequencies over the entire spectrum. 
The acoustic gap is similar in magnitude to the gaps of 
class I and spherical grains. However, there is no re- 
semblance of the triplet-doublet structure. More impor- 
tantly, Mark's decahedra do not have breathing modes. 
The small [lll]-oriented dents (see Fig. may be re- 
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FIG. 5: (a) Equivalent to Fig. but for twined icosahe- 
dral grains. Strain leads to clusters of non-degenerate fre- 
quencies, and acoustic gaps that are bigger than in single- 
crystalline grains and Mark's decahedra. (b) Mark's deca- 
hedra and multigrain's phonon frequencies up to vjj . These 
two nanograins do not have nearly-radial phonon modes. Like 
Ih's, M-Dh's show non-degenerate frequencies in this range. 
In multigrain, grain boundaries and surface disorder lead to 
non-degenerate frequencies and significantly smaller acoustic 
gap. 




Spherical 



1289 j^pj, II 








V 







i i i m~ 
1 2 3 4 5 6 7 



Multigrain 



783 




225 




n i i r 
1 2 3 4 5 6 7 



V (THz 



sponsible for these features, (iii) Multigrains only have 
non-degenerate frequencies. Like M-Dh, multigrain do 
not have breathing modes. Surface disorder and defects 
reduce the acoustic gap. (Compare Fig. |3|and0]) The 
lowest-lying modes are localized on the surface of the 
grain; however, the degree of localization is sample (dis- 
order) dependent. In addition, the more localized the 
mode the softer its frequency. The visible gap between 
the lowest allowed phonon mode and the first excited vi- 
bration in the N=459 nanogram is disorder sensitive. 



B. Vibrational density of states 



FIG. 6: Morphology and size effects on vibrational density 
of states g(v). Shaded histograms show the partial density 
of states due to surface atoms. Bin size is 0.1 THz. At low 
and high frequencies, g(v) depends significantly on geometri- 
cal details of the nanograins. Density of states at v < 2 THz 
is higher in spherical grains than in TO. For faceted and Ih 
grains, the longitudinal-phonon (LP) peaks around 4.7 THz 
have nearly the same amplitude regardless of size. In con- 
trast, LP's amplitude grows with grain size in multigrain and 
spherical grains. Grain boundaries in multigrain and twined 
Ih lead to a frequency tail beyond the maximum value in 
single-crystalline grains. No tail is seen in Mark's decahe- 
dra, whereas the LP peak is visible. [See inset in panel Ih, 
561 for g(u) of a 585 M-Dh.] Arrows indicate features arising 
from phonons strongly localized in the center of the grains. 



Figure shows the effects of morphology and size on 
the vibrational density of states (DOS) per atom g{v). In 
these results, each g{v) histogram has a 0.1 THz bin. 20 
At low and high frequencies, g(v) depends strongly on 
TO grains size. DOS in TO grains increase up to 2 THz, 
which correspond to a maximum in bulk DOS, originated 
by transversal phonons (TP). The longitudinal-phonon 
peaks around 4.7 THz have nearly the same amplitude 
for different sizes. Other class I grains (c.f. Sec. ITT}) 
show a similar behavior. DOS extends to slightly higher 
frequencies for smaller grains. Icosahedra and Mark's 
decahedra. — In icosahedra, strain and twinning remark- 
ably smear the LP peak in the DOS. Its intensity is lower 
than the TP peak about 2 THz, and remains roughly 
unchanged with increasing size. In addition, there is a 



frequency tail beyond the highest frequency in TO and 
spherical beyond 6 THz, the DOS show small features (ar- 
rows in Fig. |SJ Ih) that originate in a four phonon modes 
strongly localized at the center of the grains. These 
phonon modes consist of displacement of a few (5-40) 
atoms along the grain-boundary interfaces that meet at 
the center of grain. Removal of the center atom of the 
icosahedra results in the disappearance of these features 
in the density of states. Mark's decahedra show a DOS 
significantly different from icosahedra; see the inset in 
Fig. |SJ panel Ih, 561. Although M-Dh's are twined and 
strained, an LP peak is visible and its intensity remains 
nearly unchanged with size (not shown). More impor- 
tant, Mark's decahedra do not show a high-frequency tail 
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in DOS. Multigrain and spherical nanograms. — In spher- 
ical grains, protrusions increase g(v) below 2 and above 
5 THz relative to faceted grains. Although less promi- 
nently that in multigrains, the amplitude of LP peaks 
also increases and narrows with size. In multigrain, (i) 
the amplitude of the LP peaks grows and narrows with 
size. Remarkably, Pasquini et al 2 ^ observed a similar 
build-up of longitudinal phonons of nanocrystalline a- 
Feml (ii) Similarly to Ih, a frequency tail appears at high 
frequencies. (Arrows also point to features arising from 
inner atoms.) 

Finally, let us mention that Derlet and SwygenhovenSi 
have recently shown that the features at higher frequen- 
cies in the DOS of nanocomposites originate from vibra- 
tions strongly localized in grain boundaries. The fre- 
quency tail that appears in g(y) for Ih and multigrain 
is an intrinsic feature of nanoscopic metals with grain 
boundaries. This remarkable finding reveals the strong 
connection between extended nanocrystalline materials 
and finite size nanoscopic grains. 



C. Phonon localization 

The participation ratio N e ft(y) (see Appendix) indi- 
cates the effective number of atoms participating in a 
phonon with frequency v. Thereby, N e ff{v) measures 
the degree of phonon localization. In a sample with N 
atoms, N e ff/N ~ 1/2 for extended, wave-like phonon 
modest Note that when we identify a phonon mode as 
being localized [N e ff/N <C 1/2] we do not imply neces- 
sarily that the amplitude of the atomic vibrations decays 
away exponentially from a particular atomic site, as it 
is the case for Anderson-localized vibrational states^ 
In addition, while Anderson localization takes place in 
disordered systems, the term localization in the context 
of our results is broader as grains have a surface and, 
therefore, phonon localization — i.e., phonon modes with 
a small number of atoms involved in the vibration — takes 
place even in crystalline grains. 

Figure [7| shows N e ff{y) scaled by N, and the reduced 
vibrational density of states g(v)/v 2 for N=459 TO + , 
spherical and multigrain. The central results can be 
summarized as follows, (a) General features of N e ff 
are representative for all nanograins: (i) High-frequency 
phonons become more localized in both single-crystalline 
and multigrain specimens. In the latter, grain bound- 
aries lead to a higher degree of localization. This is 
similar to what occurs in nanocrystalline materials^ a 
result that shows that nanograins capture some of the 
physics of their massive, nanocrystalline counterparts. 26 
(ii) At low frequencies, the more localized the vibrations 
the higher the enhancement of DOS. In particular, spher- 
ical grains show a peak (arrow, Fig. 0) around 0.6 THz 
due to localized phonon modes at the surface protru- 
sions. In turn, multigrain shows a similar enhancement 
due to phonons localized around vertex vacancies and in 
atoms in the vicinity of grain boundaries. In particular, 
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FIG. 7: Morphology effects on the reduced density of states 
g{u)/v 2 and the effective number of atoms N e ff(v) in a 
phonon of frequency v. For wave-like modes, N e ff/N = 1/2. 
Only multigrain shows a v 2 distribution from 1 to slightly over 
2 THz. In all grains, edges and vertices significantly enhance 
g(v)/v 2 below ITHz. An arrow indicates modes heavily lo- 
calized on the protrusions of the spherical grain, which coin- 
cide with the peak around 0.5 THz. In the multigrain, grain 
boundaries and irregular surfaces are responsible for the re- 
duced fluctuations shown in N e f f /N. 



two low-frequency modes at 0.1 THz and 0.3 THz that 
involve ~ 5 and 50 atoms, respectively, are responsible 
for the low-frequency peak in g(y)/v 2 . Faceted single- 
crystal and twined samples (not shown) show a lesser 
enhancement, arising only from the edges and vertices, 
(iii) N e ff significantly fluctuates around N/2 for single- 
crystal samples, twined Ih and M-Dh; whereas in multi- 
grain, irregular shape, surface disorder, and defects re- 
duce these fluctuations. 

(b) Two prominent features appear in g(v)/v 2 : (i) 
The nearly quadratic behavior of g(y) from 1 — 2.2THz 
for multigrain. The frequency extent of this Debye- 
like \g{v) ~ v 2 ] region depends on the nature of dis- 
order and defects — multiple grain boundaries, stacking 
faults, and vertex vacancies. Such Debye dependence 
has been observed in iron nanocomposites by Pasquini 
and coworkers^ The Debye constant derived from the 
quadratic interval is nearly twice as big as the bulk value 
for Ag. Remarkably, such enhancement of the Debye 
constant has also been observed by Fultz and coworkers 
in nanocrystalline iron. 27 It should be noted that while 
we have found Debye-like behavior in a few nanograins 
of different sizes, the statistics we have is limited due to 
the intrinsic difficulty in generating competitive (Sec. |nj 
multigrain samples using long molecular-dynamics sim- 
ulations, (ii) The dramatic correlation between N e ff(y) 
and g(is)/is 2 , which shows that the degree of phonon lo- 
calization in nanoscopic metal grains can be derived from 
a measurement of their vibrational density of states, after 
dividing by v 2 . 

Finally, Figs. [Eta) and (b) show, respectively, displace- 
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a) Spherical b) Multigrain 




FIG. 8: Displacement vectors for two selected low-frequency, 
localized phonon modes in N=459 spherical (a) and multi- 
grain (b) samples. For each mode, we also show the frequency 
(u) and N e ff (Sec. 1111 CI and App. [XJ. In (a), the atoms 
with significant displacement-vector amplitude are those lo- 
cated on the 6-atom protrusions (indicated with solid arrows) 
that decorate the [111] facets of the underlying TO grain [Fig. 
Uc)]. Single-atom protrusions that decorate the [001] facets 
have negligible displacement-vector amplitude (open arrow). 
In (b), the phonon mode involves a very low number of atoms 
(N e ff = 6) that are around a vertex vacancy (indicated with 
an arrow) and grain boundaries. 

merit vectors for two low- frequency, localized (N e f / /N <C 
1/2) phonon modes in N=459 spherical and multigrain 
samples. These phonon modes are among those respon- 
sible for the low-frequency peaks in g(v)/v 2 (Fig. |7J). 
We find the following, (i) In the spherical grain, the 
phonon mode is one of a three-fold degenerate set with 
frequency v = 0.68 THz and N eff /N = 0.098. This 
mode is mostly localized in the 6-atom protrusions that 
decorate the hexagonal, [111] facets of the underlying TO 
grain (Fig. ^ an d Sec. [Hj. The amplitude of the displace- 
ment vectors of the atoms in the protrusions is different 
depending on which facet the protrusion is located, as in- 
dicated by solid, thick- and thin-stem arrows in Fig. Ufa) , 
whereas displacement vectors of atoms in diametrically 
opposed protrusions have the same amplitude but differ- 
ent phase. Note that in this phonon mode [Fig. Eta)] 
the single-atom protrusions located on top of [001] facets 
have negligible displacement-vector amplitude (open ar- 
rows), (ii) The phonon mode in the multigrain sample 
[Fig. EJb)] is the lowest allowed phonon. Its frequency 
is v = 0.11 THz and N eff /N = 0.013. The displace- 
ment vectors with significant amplitude locate around a 
vertex vacancy, indicated with an arrow, and along grain 
boundaries. 



IV. SUMMARY 

To understand in detail the vibrational properties of 
nanoscopic metal grains we performed atomistic simula- 
tions of the phonon spectra — frequencies and displace- 
ment vectors — of silver grains with realistic morpholo- 



gies. Our results show that morphology introduces a high 
degree of complexity into the phonon spectra, total and 
partial vibrational density of states, and phonon localiza- 
tion. The most prominent features are the following, (a) 
Low-energy, single-crystalline grains present nearly-pure 
torsional and radial phonon modes. The grains' sharp 
edges and atomic lattice split off the acoustic gap quin- 
tet predicted by ET into a doublet and triplet, with a 
magnitude that depends on the relative number of atoms 
at the boundary of different facets. When compared to 
faceted grains of the same size, high-energy, spherical 
models that present regular protrusions on the surface 
have a smaller acoustic gap and a higher total DOS at 
frequencies v < 2 THz. (b) Twined icosahedra also have 
a breathing mode. Strain in these grains leads to a higher 
acoustic gap, and a high-frequency tail in the DOS that 
originates from core atoms. This tail extends beyond the 
highest frequency in single-crystalline grains. Remark- 
ably, neither twined and strained Mark's decahedra have 
a breathing mode nor do they exhibit a high-frequency 
tail on the DOS. (c) Nanograins with grain boundaries 
and surface disorder do not have degenerate frequen- 
cies and the acoustic gap is significantly reduced. These 
nanograins are the only ones that exhibit low-frequency 
v 2 DOS in the interval l-2THz. The extent of this region 
depends on the nature of disorder. 

These predictions, while illustrated here for silver, 
should be valid for other nanoscopic metal grains. Our 
simulations show the complexity of the vibrational prop- 
erties of metal nanograins, and point out similari- 
ties between the vibrational density of states of metal 
nanograins with grain boundaries and surface disorder 
and that of massive nanocrystalline samples. Finally, 
this work is a necessary step toward understanding the 
size-dependent electron-phonon coupling in nanoscopic 
metals. 



Acknowledgments 

Funds from U. S. DOE grant DE-F G02-99ER45795 
and computational resources from the National Center of 
Supercomputer Applications (NCSA) made this research 
possible. National Science Foundation grant NSF-DMR- 
03-25939-ITR supports the Materials Computation Cen- 
ter. The authors thank Professor Jian-Min Zuo (UIUC) 
for fruitful discussions. This research started at the Ma- 
terials Computation Center (MCC). One of the authors 
(G.A.N.) acknowledges its hospitality, and the opportu- 
nity to participate in the MCC Visitor Program. 



APPENDIX A: METHODS 

Nanograin's structure and energetics. The literature 
contains extensive studies regarding the energetics of 
metal clusters. In this work, we consider known ge- 
ometries and perform total energy relaxations in Ag 
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nanograins. We also perform extensive simulated anneal- 
ing to get complex, partially disordered morphologies. 28 
To calculate the total energy S we employ Daw and 
Baskes' embeded atom method 29 (EAM). We use Voter 
and Chen's parametrization of Ag potential. 

Vibrational properties. Once we generate stable 
nanograin morphologies we calculate phonon frequen- 
cies v and displacement vectors {u*} — where i labels 
an atom located at equilibrium position R l in the 
grain — by direct diagonalization of the dynamical ma- 
trix D l ^o = d' 2 £/dR l a dR 3 j3 . R l a is a cartesian component 

(a, (3 = x,y, z) of R i . Note that we consider free-standing 
nanograins, hence, the first six phonons correspond to 
rigid-body translations and rotations along and about 



the principal axis, respectively. These frequencies do not 
enter in our analysis. 

The vibrational density of states (DOS) g{y) = 
Ef=i <?», where 9j ( V )Av = (1/N) £„ n{v)^M\Ul\ is 
the partial (local) vibrational density of states at atomic 
site j; ul is the displacement vector of atom j; n{y) is 
the number of phonons between v and v + Az/; and M is 
the mass of a Ag atom. 

The equality iV e //(z/) = M 2 /Af 2 defines the participa- 
tion ratio^Si^iiS Here, M p = JV [ei(j/)] p are momenta 
of the mean kinetic energy £i{v) of atom i. By using the 
calculated displacement vectors of atom i in a phonon 
with frequency v one readily calculates £i{v). 
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